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\d ■ ABSTRACT 

O' 

o: 

^ I We are conducting a large observing program with the Spitzer Space Telescope to 

pi ^' determine the mid-to-far infrared spectral energy distributions of a well-defined sample 

■<^ ■ of 87 nearby, 12 //m-selected Seyfert galaxies. In this paper we present the results of 

Q ' IRS low-resolution spectroscopy of a statistically representative subsample of 51 of the 

^~^ . galaxies (59%), with an analysis of the continuum shapes and a comparison of the Seyfert 

y—i . types. We find that the spectra clearly divide into groups based on their continuum 

^ ', shapes and spectral features. The largest group (47% of the sample of 51) shows very 

fSJ ! red continuum suggestive of cool dust and strong emission features attributed to PAHs. 

^ I Sixteen objects (31%) have a power-law continuum with spectral indices 05-20 //m =-2.3 

Q ' - -0.9 that flatten to 020-35 /xm =-1.1 - 0.0 at ~20 ^m. Clear silicate emission features 

^^ ■ at 10 and 18 //m are found in two of these objects (Mrk 6 and Mrk 335). A further 16% 

■^ ■ of the sample show power-law continua with unchanging slopes of 05-35 um =-1.7 - -1.1. 

o . Two objects are dominated by a broad silicate absorption feature. One object in the 

JL i sample shows an unusual spectrum dominated by emission features, that is unlike any of 

f~i [ the other spectra. Some spectral features are clearly related to a starburst contribution 

~^ I to the IR spectrum, while the mechanisms producing observed power-law continuum 

KjJ ' shapes, attributed to an AGN component, may be dust or non-thermal emission. The 

j^ ■ infrared spectral types appear to be related to the Seyfert types. Principal component 

^ ■ analysis results suggest that the relative contribution of starburst emission may be the 

. . . . dominant cause of variance in the observed spectra. The derived starburst component 

of each spectrum, however, contributes <40% of the total flux density. We compare 

the IR emission with the optically thin radio emission associated with the AGN and 

flnd that Sy I's have higher ratios of IR/radio emission than Sy 2's, as predicted by the 
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unified model if the torus is optically thick in the mid-IR. However, smooth-density torus 
models predict a much larger difference between type I's and 2's than the factor of 2 
difference observed in our sample; the observed factor of ~2 difference between the type 
I's and 2's in their IR/radio ratios above 15 /xm requires the standard smooth-density 
torus models to be optically thin at these wavelengths. However, the resulting low 
torus opacity requires that the high observed columns detected in X-ray absorption be 
produced in gas with very low dust to gas ratio (perhaps within the dust sublimation 
region). On the other hand, our observations may be consistent with clumpy torus 
models containing a steep radial distribution of optically thick dense clumps. The 
selection of our sample at 12 fira, where the torus may be optically thick, implies that 
there may be orientation-dependent biases in the sample, however we do not find that 
the sample is biased towards Sy 2's with more luminous central engines as would be 
expected. We find that the Sy 2's typically show stronger starburst contributions than 
the Sy I's in the sample, contrary to what is expected based on the unified scheme for 
AGN. This may be due to the selection effect that only those Seyfert 2's with strong 
starburst contributions had high enough integrated 12 //m fiux densities to fall above 
the flux limit of the sample. 

Subject headings: galaxies: Seyfert — galaxies: spiral — infrared: galaxies 



1. Introduction 

Active galactic nuclei (AGNs) are powered by accretion onto a massive black hole (see reviews 
by Rees 1984 and Shields 1999). The fuel itself is probably ISM from the host galaxy that is driven 
to the center by bar-induced torques (Jogee et al. 2005), galaxy interactions, or galaxy mergers 
(Osterbrock 1993 and references therein). Starbursts may be an inevitable consequence of the gas 
infall (Rees 1984; Norman &: Scoville 1988), and the observed relation between black hole mass 
and host bulge mass in AGNs suggests that nuclear activity is closely connected to star formation 
(Merritt & Ferrarese 2001; Gebhardt et al. 2000). Nuclear activity, though probably short-lived 
compared with the lifetime of the host (~5 - 20% percent of galaxies host AGN; Maiolino & 
Rieke 1995; Ho et al. 1997; Kauffmann et al. 2003; Hao et al. 2005a), may be an integral part of 
the evolution of all galaxies. Understanding the inner workings of AGNs and the AGN/starburst 
connection is therefore important for understanding galaxy evolution. 

Seyfert galaxies are the nearest and brightest AGNs and so are well-suited for studying the 
connection between nuclear activity and star formation. Obscuring dust however hampers studies 
of the optical - soft X-ray emission produced by young stars and the accretion disk. Emission 
at infrared wavelengths does not suffer such large extinction, and, further, the dust that absorbs 
the shorter wavelength emission reradiates in the IR and correspondingly produces a substantial 
fraction of the bolometric flux of the object. Dissecting the detailed IR spectral energy distribution 



(SED) of Seyfert galaxies can reveal the properties of the dust in the nuclear region as well as place 
constraints on the optical - soft X-ray spectrum (i.e., starburst and AGN spectrum) that heats the 
dust. 

A key factor in understanding Seyfert galaxies is determining the geometry of the nuclear 
obscuring material. Unification schemes postulate that the observed differences between Seyfert 
types 1 and 2, such as optical emission line widths and X-ray spectral slopes, are due to orientation- 
dependent obscuration (see reviews by Urry & Padovani 1995; Antonucci 1993). Alternatively, 
there may be fundamental differences between the central engines of Sy I's and 2's. The obscuring 
medium has been modeled as a dusty torus which, viewed edge-on (Sy 2's), obscures the nucleus 
but, viewed face-on (Sy I's), allows the nucleus to be seen. The inner scale of the torus is thought 
to be of order 1 pc (e.g., Krolik & Begelman 1988; Gallimore et al. 1999; Ulvestad et al. 1999; Jaffe 
et al. 2004; Rottgering et al. 2005), but other properties are poorly constrained owing to angular 
resolution limitations. However, the geometry and equatorial opacity of the torus should have an 
observable impact on the shape of the infrared SED (Pier & Krolik 1992; Granato & Danese 1994; 
Ivezic & Elitzur 1997; Nenkova et al. 2002). Thus AGN unification can be tested by comparing the 
SEDs of type 1 and 2 Seyferts with model predictions. 

A third possibility is that the Seyfert types are connected by evolution: an active nucleus 
evolves through a starburst (type 2) phase during which it is heavily obscured, and into an AGN 
phase (type 1) when the obscuring dust has been evacuated from the nucleus by starburst-driven 
winds and AGN outflows (Sanders et al. 1988; Storchi-Bergmann et al. 2001). Supporting this 
picture is evidence that Seyfert 2's have significantly more star formation than Seyfert I's. Photo- 
metric studies have found that Sy 2's have fainter nuclear infrared emission than Sy I's and SEDs 
consistent with starburst-dominated IR emission, while Sy I's have SEDs consistent with pure AGN 
+ quiescent galaxy emission (Edelson, Malkan, & Rieke 1987; Maiolino et al. 1995). However, as 
orientation is also likely to play a role in the differences between Seyfert Is and 2s, it is important 
to investigate the relative contributions of star formation and the active nucleus to the infrared 
emission in order to test both AGN unification and evolutionary schemes. 

The Spitzer space telescope represents a significant step forward in both sensitivity and spatial 
resolution compared with previous thermal infrared missions. We are conducting a large observing 
program with the Spitzer space telescope to determine the mid-to-far infrared (MFIR) spectral 
energy distributions (SEDs) of 87 Seyfert galaxies (PI: J. Gallimore, PID: 3269). The majority of 
the IRS spectra have been observed (51 objects, 59%), allowing the first statistical analysis of the 
detailed mid-infrared (MIR) SEDs of a relatively unbiased sample of Seyfert galaxies. In this paper 
we present an overview of the continuum properties of the sample, and relate these to the optical 
spectral types. An atlas of spectra and detailed modeling of individual spectra will be presented in 
a future paper when all data are available. The purpose of this paper is to present the results of our 
initial analysis of the spectra in order for them to be useful to the AGN community. In §2 we briefly 
describe the sample properties. The observations and the data reduction are presented in §3. §4 
presents typical spectra and an analysis of the continuum shapes. In §5 we discuss the implications 



of our results for the unification of Seyfert galaxies and our understanding of the nature of the 
torus. 



2. The Seyfert sample and its properties 

The sample of Seyfert galaxies used for this study comprises all Seyfert galaxies from the 
extended 12 fira sample of Rush et al. (1993) that have cz < 10000 kms~^. The potential advantage 
of the 12 //m sample is its perceived resistance to wavelength-dependent selection effects (Spinoglio 
&: Malkan 1989). The redshift criterion was included to ensure that the smallest possible physical 
size was probed by the fixed Spitzer apertures in order to better exclude host galaxy emission. 
The resulting sample contains 87 sources. Three of the sample galaxies, NGC 1097, NGC 1566, 
and NGC 5033, are being extensively observed as part of the SINGS Legacy project, so we are 
obtaining data for 84 galaxies. The data include IRAC imaging, IRS spectral mapping, and MIPS 
SED spectra. The present paper discusses the 51 objects in the sample for which 5-35 /im IRS 
spectra are currently available. 

The spectral classifications were obtained from the Veron-Cetty &: Veron (2003) catalog, and 
thus have not been made using homogeneous datasets of consistent quality. In several cases, the 
classification is ambiguous, and we have adopted revised classifications from other sources (Tran 
2003; Ho et al. 1997; Veilleux et al. 1995; Phillips 1979). Seven objects in the sample have been 
reclassified as non-Seyferts (LINERs or star-forming galaxies): MCG -hO - 29 - 23, NGC 4922, 
NGC 5953, NGC 7130, NGC 7496, NGC 7590, and UGC 11680. These objects have been excluded 
from analyses comparing Sy I's and Sy 2's, and we expect the results will not be sensitive to the 
classifications of these objects. For the purposes of comparing the Seyfert types we consider Seyfert 
types l.n (i.e. types 1.0, 1.2, 1.5, 1.8, and 1.9 as well as narrow line Seyfert I's) to be Sy I's, and 
objects with no broad permitted lines or with hidden broad line regions (HBLRs) to be Sy 2's. 

Figures 1 and 2 show the distributions of the IRAS 12 /im fiux densities and redshifts of the 
sample objects. Data for the whole sample are shown, with the objects observed to date with 
Spitzer (~ 50% of the sample) indicated by hatching for comparison. Sy I's and Sy 2's are shown 
separately. Two-sample statistical tests were performed to determine whether the distributions of 
flux density and redshift are similar for the Seyfert types and if the objects observed to date are 
representative of the 12 ^um sample as a whole. The probabilities of the various datasets being 
drawn from the same parent population were calculated using the ASURV package (Lavalley et al. 
1992), which enables the analysis of data in which limits are present (Feigelson &: Nelson 1985). 
The results of Gehan and logrank tests are shown in Table 1. We find no significant differences 
between Sy I's and Sy 2's in either the IRAS 12 ^um flux density or redshift distributions. We also 
find no significant differences between the observed subsample and the whole sample, indicating 
that the current dataset is broadly representative of the complete sample. 

The optically thin nuclear radio emission may be considered an isotropic property of Seyfert 
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Fig. 1. — Distribution of the IRAS 12 yum fluxes for the Seyfert I's (top) and the Seyfert 2's 
(bottom). Hashed regions indicate the objects with IRS data in hand. Two-sample statistical tests 
show no significant diff'erence between the Sy I's and Sy 2's or between the observed objects and 
the whole sample (Table 1). 
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Fig. 2. — Distribution of redshifts for the 12 micron sample. Seyfert I's (top) and Seyfert 2's 
(bottom) are shown separately. Hashed regions indicate the objects with IRS data in hand. Two- 
sample statistical tests show no significant difference between the Sy I's and Sy 2's or between the 
observed objects and the whole sample (Table 1). 



galaxies, that provides an orientation-independent measurement of the energy output of the active 
nucleus (Xu et al. 1999; Edelson 1987; Giuricin et al. 1990; Kukula et al. 1995; Thean et al. 2001). 
Radio observations of the 12 fiui Seyfert sample have been obtained by Thean et al. (2000), using 
the VLA A array at 8.4 GHz. These high spatial resolution (0.25"), high frequency observations 
isolate the nuclear radio emission (< 3.5" in extent) associated with the AGN from radio emission 
associated with star formation in the host galaxy. Seyfert nuclei are relatively weak radio emitters, 
so the larger scale radio emission associated with star formation can dominate over the nuclear 
radio emission at lower spatial resolutions. Figure 3 shows the distributions of nuclear radio flux 
densities for the whole sample and for the objects for which we have IRS data in hand. Sy I's 
and Sy 2's are shown separately. Limits are shown for 4 objects not detected by Thean et al. 
(2000). Radio flux densities for 6 objects not observed by Thean et al. (2000) were obtained from 
NED. These were generally lower resolution and lower frequency, and so upper limits to the flux 
density at 8.4 GHz were derived assuming Si, oc v^^'"^ . Wc find no significant difference between 
the distributions of nuclear radio flux density for the type 1 and 2 Seyferts. Results of statistical 
comparisons of sample properties are summarized in Table 1. 



3. Observations and Spectral Extraction 

The sample galaxies were observed by Spitzer IRS (Houck et al. 2004) in mapping mode, using 
the Short-Low and Long-Low modules. The modules cover a wavelength range of --^5 - 35 /xm, 
with a resolving power of 64 - 128. The raw data were processed by the Spitzer pipeline version 
11.0. The nuclear spectra were extracted from the Basic Calibrated Data (BCDs) using SMART 
(Higdon et al. 2004). Preliminary spectral maps show that the surface brightness is dominated by 
a central point source at these wavelengths, so the default point source aperture was used for the 
extractions. Only three of the sources (NGC 1365, NGC 5005, &: NGC 5953) show a noticeable drop 
in flux between the overlapping LL (long-low) and SL (short-low) module spectra, which results 
from extended flux contributing to the wider LL module slit (SL: 3.5", LL: 10.5"). This supports 
our conclusion that the objects are dominated by the central point source, and so default point 
source extractions are appropriate. Flux densities were calibrated based on flux calibration tables 
provided by the Spitzer Science Center. The spectra are shown in Figure 4, grouped according to 
their shape (see §4.2). The full spectral maps will be presented in a future paper. 

Following a suggestion from the referee, comparison of the flux densities of the nuclear IRS 
spectra with IRAC 5.8 and 8.0 ^m nuclear photometry (Gallimore et al. 2005, in preparation) later 
revealed that the spectra underestimated the source flux densities. We attribute this to slit losses 
that vary as a function of wavelength due to the increasing size of the point spread function with 
increasing wavelength. As these spectra were obtained in mapping mode, no attempt was made to 
center the nuclear point source in the slit during the observations. Thus the nuclear (central) slit 
position is not necessarily well centered on the source and the spectra suffer wavelength-dependent 
slit losses, up to 20%. This is only a problem for strongly point-source dominated sources. We have 
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Fig. 3. — Distribution of nuclear radio flux densities for the 12 micron sample. Seyfert I's (top) 
and Seyfert 2's (bottom) are shown separately. Hashed regions indicate the objects with IRS data 
in hand and arrows indicate upper limits. The results of two-sample tests comparing distributions 
are shown in Table 1. 



attempted to account for this by producing a spectral image cube from the Basic Cahbrated Data 
FITS files of all the slit positions (after background subtraction using the off orders) and iteratively 
determining an interpolated nuclear image based on the fractional shift perpendicular to the slit 
that produced the most total flux in the interpolated image. A simplex maximization technique 
was used to find the optimum interpolated position, which defined the plane containing the spectral 
image of the point source. We refer to the spectra extracted from the spectral image of the point 
source as "pointing-corrected" nuclear spectra. Flux densities derived from the resulting pointing- 
corrected spectra agree well with our IRAC photometry, within 20% for the vast majority of objects 
and 30-45% for a few objects. The IRS-derived flux densities both overestimate and underestimate 
the IRAC 5.8 and 8.0 //m point source flux densities, indicating that the flux calibration of our IRS 
spectra is not systematically underestimating the true flux densities. Excluding the comparison 
with IRAS (§4.1), in the remainder of the paper the analysis was performed on the original (not 
pointing-corrected) spectra. Despite the fact that the flux errors are wavelength-dependent, we 
find that the overall shapes of the spectra are not changed significantly by the pointing correction. 
We do not expect that the flux errors will have a significant impact on the results. 



Table 1. Properties of the sample: results of two-sample tests. 
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Note. — Generalized Wilcoxon and Logrank two-sample tests were used to determine the probability that the values of the variable (column 
[1]) in groups 1 and 2 (columns [2] and [3]) were drawn from the same parent distribution. The number of objects in the groups are: all Sy I's: 
37; observed Sy I's: 22; all Sy 2's: 38; observed Sy 2's: 22. Columns [4]-[6] list the parameters for group 1: the number of limits in the data, 
the mean value and the median value, derived using the Kaplan-Meier (K-M) estimator. * Indicates the mean is biased because the lower- or 
uppermost limit was changed to a detection to compute the K-M distribution. Columns [7]-[9] list the same parameters for group 2. Large values 
of p (columns [10] and [11]) indicate there is no significant difference between the distributions for the two groups. Column [12] refers to the 
Figure number showing the distribution in the variables. 
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Fig. 4. — (a) Spectra dominated by red continuum and PAH features. Each spectrum is labeled 
with the object name and the flux density offset applied. 
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Fig. 4. — (b) The left panel shows spectra that may be described by a broken power law continuum. 
The right panel shows spectra that can be fitted with a continuous power law, the two spectra 
showing silicate absorption at 10 //m, and the ungrouped spectrum (bottom). Each spectrum is 
labeled with the object name and the flux density offset applied. 
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4. Results 

4.1. Comparison with IRAS 

The IRS apertures are much smaller than those of IRAS, and so we can evaluate the fraction 
of extended flux not included in the IRS aperture (assuming insignificant mid-infrared variation in 
the ~ 24 years between observations). The IRS spectra span the wavelength range 5-35 /xm, so 
12 and 25 /^m nuclear flux densities corresponding to the IRAS bandpasses can be derived from the 
nuclear spectra. For the purposes of this comparison, the nuclear flux densities were determined 
using the pointing-corrected IRS spectra and the IRAS bandpass spectral response functions. The 
distributions of nuclear 12 and 25 ;um flux densities for Sy I's and Sy 2's are shown in Figure 
5. There is a marginal probability that the distributions for Sy I's and Sy 2's are different at 
12 una, with the Sy 2's fainter than the Sy I's (7% probability that the difference in the observed 
distributions ocurrs due to chance; Table 2), but the distributions are statistically similar at 25 ^m. 

We calculate the flux deficit, the fraction of flux "missing" in the IRS aperture relative to 
IRAS, Fdef = {FiRAS — FiRs)/FiRAS- Figures 6 and 7 show that the sample includes some 
objects for which all the IRAS flux is concentrated within the IRS low-resolution slits (widths 3.5 
- 10.5"), and other objects for which most of the flux falls outside the slit. For most objects, 
the majority of the IRAS flux is seen by IRS, indicating the surface brightness distributions of 
the objects are compact. For a few objects, the IRS flux density is higher than the IRAS flux 
density, suggesting that variability may be contributing to the scatter in the distributions of flux 
deficit; however, the differences are <20%, which is within our flux density uncertainties. These 
results are consistent with the results of previous studies that compared ground-based photometry 
of Seyfert galaxies with IRAS photometry. Maiolino et al. (1995) compared 5" aperture 10 /im flux 
densities of Seyfert galaxies with the IRAS measurements and showed that the Seyferts displayed a 
range of ground/IRAS flux densities from 0.0 - 1.0, indicating that the sources ranged from highly 
compact objects dominated by the nuclear point source, to those with relatively strong extended 
emission, associated with star-formation. Gorjian et al. (2004) also compared small aperture (1.5") 
ground-based photometry of Seyfert galaxies with IRAS measurements. Their results show a similar 
distribution of flux deficits to our dataset, and these authors note that their data are relatively 
insensitive to diffuse extended emission due to the high spatial resolution of the detector, and to 
the small chopping distance used in the observations, which likely resulted in over-subtraction of 
background emission. 

We find that the average fiux deficits of Sy 2's (means of 0.71 and 0.66 at 12 and 25 ^m, 
respectively) are greater than those of Sy I's (means of 0.62 and 0.56; Figures 6 and 7), although 
statistical comparisons indicate that the differences are only marginally significant (formal proba- 
bilities that the parent distributions are the same are 6 - 11%; Table 2). The lack of a jump in 
flux density between the spectra obtained with the SL and LL modules, which have different slit 
widths, suggests that the emission occurs on scales larger than the LL module slit width (~10"). 
This implies that much of the MIR emission is on kiloparsec scales, and that this difference between 
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the Sy I's and 2's is related to relatively extended star formation rather than orientation-dependent 
obscuration in the nucleus. The difference in flux deficits between the Sy I's and Sy 2's implies that 
the Sy 2's have relatively more extended emission than the Sy I's, that is, they are less dominated 
by the nuclear point source. If both types of Seyfert have the same luminosities, then our results 
would indicate that the type 2's have more luminous extended emission than the type I's. As 
noted above, the distributions of nuclear 25 /um flux density for the Sy I's and Sy 2's are similar, 
while at 12 ^um the nuclear flux densities of type I's and 2's are marginally significantly different 
(Fig. 5 and Table 2). We note that, as the sample was selected at 12 ^um, Sy 2's with weak nuclei 
may have been excluded, resulting in the 12 /xm distributions of type I's and 2's appearing more 
similar than they in fact are; the effects of sample selection on our results are discussed further in 
§5.3. Therefore it would appear that the Sy 2's have more extended emission relative to the nuclear 
emission than the Sy I's and that the Sy 2's have, on average, more luminous extended emission 
than the type I's. 



4.2. Spectral shapes 

We find four distinct types of continuum shapes and spectral features among the Seyfert nuclear 
spectra. Figure 8 shows typical spectra in each group. An atlas of spectra will be presented when 
all data have been obtained. All spectra have been converted to the restframe. We group the 
spectra according to the following main properties: 

• Twenty-four objects (47% of the observed objects) have PAH-dominated spectra with strong 
emission features at 6.2, 7-9, 11.3, and 12.9 fim and very red continuum suggestive of cool dust; 
Mrk938 is the archetype (Fig. 8). Eight objects, including Mrk938, show clear silicate absorption 
at 10 /um, while a further eight show possible weak absorption, but the continuum level is difficult 
to determine due to the strong PAH emission. The remaining 8 objects in this group show no 
apparent silicate features. There are 5 Sy I's, 12 Sy 2's, 3 LINERs, and 4 starburst galaxies in this 
group. 

• Sixteen objects (31%) have continua that can be described by a broken power-law; they show a 
flattening in the continuum slope at ~20 //m; NGC4151 is the brightest of this class (Fig. 8). This 
flattening may be due to a warm (--^170 K) dust component, peaking at -^20 ^m, dominating the 
mid-IR emission (e.g., Weedman et al. 2005; Perez Garcia &: Rodriguez Espinosa 2001; Rodriguez 
Espinosa et al. 1996), however, a simple model of an emissivity- weighted blackbody function does 
not fit the continuum of these spectra well. Models including multiple dust components are required 
to determine if this is responsible for the continuum shape (Buchanan et al., in prep.). Two of 
these objects, Mrk6 and Mrk335, show clear silicate emission features (Figure 8b). A further 9 
objects in this group, including NGC4151 (Fig. 8), show weak excesses at 10 and 18 /um that may 
be due to silicate emission; the features are sufficiently weak that a more quantitative analysis is 
necessary to confirm the presence of silicate dust emission in these spectra. One object, NGC5347, 
shows apparent weak silicate absorption at 10 fim in its spectrum, while the remaining 4 objects 
have ambiguous spectral features. There are 9 Sy I's, and 7 Sy 2's in this group. Power-law fits to 
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Fig. 5. — Distribution of the 12 fim (left) and 25 fiuifright) nuclear fluxes measured from the IRS 
spectra. Seyfert I's (top) and Seyfert 2's (bottom) are shown separately. 
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Fig. 6. — IRAS global flux density vs IRS nuclear flux density at 12 ^um (left) and 25 nni(right). 
Seyfert I's (open circles) and Seyfert 2's (filled circles) are shown separately. The dashed lines 
indicate equal global and nuclear fluxes. 



Table 2. Properties derived from IRS spectra: results of two-sample tests. 



Variable 


Group 1 


Group 2 




Group 


1 




Group 


2 


P 


P 


Figure 








"Km 


Mean 


Median 


"Km 


Mean 


Median 


(Gehan) 


(Logrank) 




[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


IRSlog(Fi2/im) (Jy) 


obs Sy I's 


obs Sy 2's 





-0.76 


-0.78 





-0.99 


-0.92 


0.07 


0.07 


5 


IRSlog(F25/,m) (Jy) 


obs Sy I's 


obs Sy 2's 





-0.41 


-0.47 





-0.52 


-0.47 


0.78 


0.44 


5 


12 /^m flux deficit 


obs Sy I's 


obs Sy 2's 





0.62 


0.57 





0.71 


0.70 


0.06 


0.11 


7 


25 fim flux deficit 


obs Sy I's 


obs Sy 2's 





0.56 


0.57 





0.66 


0.63 


0.06 


0.07 


7 


E/v 1 component 


red, PAH 


broken P-L 





14.1 


14.1 





-7.3 


-7.3 


0.0000 


0.0000 


13 


E/v 1 component 


red, PAH 


unbroken P-L 





14.1 


14.1 





-1.1 


-2.1 


0.0001 


0.0000 


13 


E/v 1 component 


broken P-L 


unbroken P-L 





-7.3 


-7.3 





-1.1 


-2.1 


0.0003 


0.0002 


13 


E/v 1 component 


obs Sy I's 


obs Sy 2's 





-1.2 


-2.9 





6.2 


3.3 


0.04 


0.02 


14 


log(Fio.6 ^m/S'8.4 GHz) 


obs Sy I's 


obs Sy 2's 


5 


2.01* 


1.91 


2 


1.16* 


1.25 


0.002 


0.011 


16 


E/v 1 component 


high L8.4GHZ 


low L8.4GHZ 





5.6 


0.7 





4.7 


2.9 


0.88 


0.59 


17 



Note. — Generalized Wilcoxon and Logrank two-sample tests were used to determine the probability that the values of the variable (column [1]) in 
groups 1 and 2 (columns [2] and [3]) were drawn from the same parent distribution. The number of objects in the groups arc: all Sy I's: 37; observed 
Sy I's: 22; all Sy 2's: 38; observed Sy 2's: 22. Columns [4]- [6] list the parameters for group 1: the number of limits in the data, the mean value and 
the median value, derived using the Kaplan-Meier (K-M) estimator where there are limits in the sample. * Indicates the mean is biased because the 
lower- or uppermost limit was changed to a detection to compute the K-M distribution. Columns [7]- [9] list the same parameters for group 2. Large 
values of p (columns [10] and [11]) indicate there is no significant difference between the distributions for the two groups. Column [12] refers to the 
Figure number showing the distribution in the variables. 
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these spectra indicate that the spectral indices below 20 fim. range from -2.3 to -0.9, with a median 
value of a5_20 = -1-3 (where S^ oc z^"). The spectral indices above the break are -1.1 - 0.0, with 
the median 020-35 = -0-4. 

• Eight objects (16%) have power-law continuous spectra over the IRS spectral range ~5 - 35 //m. 
NGC 3516 is representative of this class (Fig. 8). These spectra show no strong dust emission or 
absorption features, though most have weak excesses at 10 and 18 ^m that may be due to silicate 
emission (Fig. 8). Six of these objects are Sy I's, and two are Sy 2's. Power-law fits to these data 
indicate the spectral indices range from -1.7 to -1.1, with a median of 05^35 = -1.2. 

• Two objects (NGC 1194 Sz F04385— 0828) show a broad absorption feature at 9.7 /um due to 
silicate dust. NGC 1194 is a Sy 1, and F04385-0828 is a Sy 2. 

• One object (NGC 7603) appears to show both PAH and silicate emission features and otherwise 
fails to meet any of the above classifications. This object is a Sy 1. 

The spectral shapes may be quantitatively characterized by the IR colors derived from the 
spectra. We calculate the spectral index between two wavelengths using 

a2-i = \og{fxJfx^)/log{Xi/X2). 

The spectral index between 20 and 30 /um quantifies the slope of the spectrum (above the -^20 /um 
'break' for broken power-law spectra), while the spectral index between 8 and 10 ^m indicates 
the presence or absence of the 9.7 ^um silicate feature and PAH features. The resulting color-color 
diagrams are given in Figure 9. One Sy 2 (NGC3079) has extreme colors 020-30 ^um = 4.7 and 
as-io/um = —12.5, due to its very red spectrum and strong PAH emission, so it is not shown 
in Fig. 9. The Sy I's (Seyfert types 1.0, 1.2 and 1.5) lie in the lower-right part of the diagram, 
occupied by unbroken and broken power-law spectra, specifically avoiding the region of the reddest 
PAH emitters. Sy 2's are found in both regions. Earlier type Seyferts (i.e. 1 - 1.5) have the bluest 
colors, while Seyfert 2's with hidden broad line regions (HBLRs; type 2h) span the range of colors. 
The fact that Sy I's consistently show the bluest infrared colors and the Sy 2's are the reddest 
sources indicates the IR spectra are changing systematically with Seyfert type. We consider the 
implications further in § 5. 
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Fig. 7. — Distribution of the flux deficit, the fraction of flux missing in the IRS spectra at 12 ^um 
(left) and 25 fim (right) compared with the IRAS flux densities. 
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been converted to restframe. 
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Fig. 9. — IR color-color diagrams derived from the IRS spectra, (a) The symbols indicate the shape 
of the spectrum: red continuum with PAH features (stars), broken power-law (triangles), unbroken 
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three largest groups; no ellipse is shown for the 2 objects with silicate absorption, (b) The numbers 
indicate the Seyfert spectral type. Numbers 1.0 - 1.9 and In represent Seyfert 1 subtypes. Symbols 
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the letter H. 
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4.3. Principal Component Analysis 

We performed a principal component analysis (PCA) on the spectra to determine the compo- 
nent spectral shapes (eigenvectors) producing the variety of spectra seen in the sample (see, e.g., 
Francis & Wills 1999; Shang & Wills 2004). Because we are presently more interested in the shape 
of the mid-infrared SED rather than narrow emission features, we median-smoothed all of the spec- 
tra to improve the continuous signal in the primary eigenvectors. We adopt a power law spectrum 
with index a = —1.2, (the average slope of the power-law spectra, where S^ oc v"-) as the AGN 
'mean' spectrum, rather than using the average spectrum derived from the data, to try to separate 
this component of the spectra from other features, although using the actual mean spectrum gives 
similar results. The spectra were converted to the restframe and normalized to unity flux density 
at 19.0 iiTQ. before performing the analysis. The PCA computes n eigenvectors, where n is the 
number of wavelength bins, but for the present data only the first few contribute significantly to 
the variance. In fact, the first eigenvector is the most dominant and contributes 91% of the variance 
in the spectra. The first 3 eigenvectors together account for 97.7% of the variance, and the first 6 
eigenvectors 99.0%. Figure 10 shows the power law spectrum and the first three eigenvectors. 

The first eigenvector shows broad emission features, a deep absorption feature at ~10 ;um, and 
red continuum above 20 /im. We interpret the first eigenvector to be dominated by a starburst 
component, as it closely resembles the spectra of starburst galaxies (Fig. 11). The emission features 
are attributable to PAH emission, while the absorption feature is due to silicate dust absorption 
accentuated by the surrounding PAH emission. The change in continuum slope at ~20 ^um in the 
PCA eigenvector 1, which is not seen in starburst spectra, probably arises from the broken power- 
law spectra in the sample, which have negative weights for this eigenvector, and suggests that the 
first eigenvector does not represent a 'pure' starburst. The second eigenvector is dominated by 
PAH features and so accounts for differing relative strengths of the PAH emission with respect to 
the silicate absorption in the spectra. The third eigenvector appears to show narrow line emission, 
in particular [O IV] A25.9 ^m, indicating that the smoothing did not completely remove these lines 
from the spectra. 

The individual spectra, or component parts thereof, can be reproduced by summing the 'mean' 
spectrum and the eigenvectors multiplied by their respective weights (e.g., Shang & Wills 2004). 
The absolute and relative contributions of each eigenvector to each spectrum can be determined 
using the associated weights. The relative contribution of the first eigenvector to each spectrum 
appears to be related to the 12 ^um flux deficit (Figure 12), supporting the conclusion that it 
predominantly indicates the starburst component of the SED. This figure further suggests that 
objects with more extended starburst activity also have more nuclear starburst activity. Figure 13 
shows that the relative contribution of eigenvector 1 to each spectrum relates strongly to the shape 
of the SED, as the different shapes have significantly different contributions from this eigenvector 
(probabilities <0.03% that the distributions are from the same parent distribution, see Table 2). 
Negative relative contributions indicate that this component is subtracted from the 'mean' spectrum 
to produce the observed spectrum. Negative eigenvector 1 components in broken and unbroken 
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Fig. 10. — The adopted 'mean AGN' power-law spectrum (top) and first three eigenvectors of the 
sample of spectra, derived through principal component analysis. 
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power-law spectra may incorporate 10 fim. silicate dust emission. The spectra with red continuum 
and PAH emission features show the largest eigenvector 1 components, consistent with these objects 
having the greatest starburst contribution. 

The relative contribution of the first eigenvector differs between Seyfert types 1 and 2, with 
Sy I's having a negative median component and Sy 2's having a positive median contribution. This 
supports the result shown in Figure 9 that the spectra showing red continua and PAH emission 
tend to be Sy 2's while the Sy I's tend to have other shapes. However, Figure 14 also shows 
that the Seyfert types have overlapping distributions of the first eigenvector contribution. The 
differences in the first eigenvector contribution are less significant for Seyfert type than for SED 
shape. Consistent with this is the lack of separation of Seyfert types in Figure 12. The statistical 
probabilities determined from two-sample tests are listed in Table 2. We are currently performing 
radiative transfer modeling of the dust emission in order to determine the physical mechanisms 
producing the different IR shapes and compare the Seyfert types in more detail (Buchanan et al., 
in prep). 



24 



100.00 




Wavelength (micron) 

Fig. 11. — The first eigenvector of the PC A which we interpret as dominated by a starburst 
component but also including contribution from other features (thin line), and, for comparison, 
the ISO spectrum of starburst galaxy M82 (thick line). The ISO spectrum was obtained from 
http://isc.astro.cornell.edu/ sloan/library/swsatlas/atlas.html. The eigenvector clearly shows the 
PAH features and silicate absorption seen in starburst galaxies. The shape of the eigenvector at 
the longer wavelength differs slightly from a starburst shape. It is possible that this shape in the 
eigenvector produces the break in slope in the NGC 4151-like objects. 
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Fig. 12. — The flux deficit at 12 //m (see §4.1) compared with the relative contribution of the first 
eigenvector derived from the PCA. The symbols indicate the Seyfert type: Sy 1 (open circles), Sy 2 
(filled circles), LINER (open square), and starburst galaxy (open star). 
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4.4. IR/Radio flux ratios 

The unified scheme for active galaxies attributes the differences between Sy I's and Sy 2's to 
the presence of an opticahy thick dusty torus surrounding the nucleus (see reviews by Antonucci 
1993; Urry & Padovani 1995). Models of such a torus predict that, if the dust is optically thick at 
mid-IR wavelengths, Sy I's (face-on) will be stronger mid-IR emitters than Sy 2's (edge-on) (Pier 
&; Krolik 1992; Granato & Danese 1994). Heckman (1995) tested this prediction by comparing the 
mid-IR emission of a heterogeneous sample of Seyfert galaxies. He derived the ratios of 10.6 /um 
to non-thermal radio flux densities, to normalize the IR emission from each object with respect to 
its intrinsic AGN brightness, and compared the IR/radio flux density ratios of the type 1 and type 
2 objects. As noted in §2, the optically thin nuclear radio emission associated with the AGN is 
isotropic and considered to be representative of the intrinsic AGN power. Heckman (1995) found 
that Sy I's are indeed stronger 10.6 /xm emitters than Sy 2's, relative to their isotropic "intrinsic" 
AGN brightness, by a factor ~4. This result extended the work of Maiolino et al. (1995), who found 
that the Sy I's have brighter absolute nuclear IR emission than Sy 2's, showing that the difference 
is due to anisotropy in the mid-IR emission. We derive the nuclear IR to radio flux density ratios 
of our homogeneous sample of Seyfert galaxies, using our IRS spectra and the nuclear radio flux 
densities (Fig. 3). Figure 15 shows the distributions of the ratio 6*10.6 /im/5'8.4GHz for the Sy I's 
and Sy 2's. We find the median of the Sy I's to be higher than that of the Sy 2's, by a factor ~5 
(see Table 2). 

We repeated this analysis at intervals of 0.5 - 1.0 //m along the IRS spectral range, in order 
to investigate this difference between Sy I's and Sy 2's across the mid-IR. The median ratio of 
log(5iR/58.4GHz) is consistently higher for Sy I's than for Sy 2's at all mid-IR wavelengths (Figure 
16), though the statistical significance of the difference decreases at longer wavelengths. The 
difference between Sy I's and 2's is particularly apparent around 10 fun. This may be due to 
strong PAH emission and/or silicate absorption in the Sy 2's that is absent in the Sy I's. Figure 16 
shows the ratio of 5'iPi,/58.4GHz of Sy I's to that of Sy 2's. The sharp increase in this ratio towards 
shorter wavelengths might indicate the presence of a hot dust component, peaking ~2 /xm, in the 
Sy I's that is absent in the Sy 2's. A higher ratio of IR/radio emission in Sy I's is consistent with 
AGN unification, if the obscuring torus is optically thick in the IR, however a clear separation of the 
starburst and AGN components in the spectra, from detailed modeling, is necessary before further 
interpretation of Figure 16 will be possible (see §5.2). In particular, it is unclear if the 10 /im peak 
arises in the torus (Sy 2's being edge-on) or in the starburst component (preferentially stronger in 
Sy 2's). 
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Fig. 15. — Histograms of the 10.6 /iin to 8.4 GHz flux density ratio for the Sy I's (top) and Sy 2's 
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Discussion 



We now consider the implications of our results for our understanding of the starburst/AGN 
connection, the AGN unified model, and selection effects in AGN samples. 



5.1. The relative IR contributions of starburst and AGN 

Spitzer spectra of the sample reveal distinct continuum shapes and dust features. The similarity 
in the features between the red objects like MRK 938 (see Fig. 8) and the IR spectra of starburst 
galaxies such as M 82 (see Fig. 11) strongly suggests these objects are dominated by the starburst 
contribution to the dust heating at mid-infrared wavelengths. We find no difference in the redshifts 
of these objects compared with the other IR spectral types, indicating that the relative strength 
of the starburst component is not an aperture effect. We do find evidence, however, that the 
strength of the nuclear starburst, as measured by the contribution of the first eigenvector in the 
PGA, is related to the extended star formation, measured by the 12 fim flux deficit (Figure 12). 
In the starburst-dominated objects the dust is cool and so peaks longward of 30 /xm, and the 
PAH features are strong. In contrast, objects with broken and unbroken power-law spectra are 
probably dominated by the AGN contribution to the spectrum. Previous studies suggest that this 
AGN contribution is a warm dust component (e.g., Rieke 1978; McAlary & Rieke 1988; Barvainis 
1987; Pier & Krolik 1993; Polletta et al. 2000). It is possible that the power-law spectra are 
produced by dust with a range of temperatures (Panagia &: Felli 1975), and the apparent break 
in the broken power-law spectra is produced by a dominant dust component peaking at ~20 fim. 
MFIR photometric SEDs of several Seyferts, based on ISO data, are consistent with such a dust 
component (Rodriguez Espinosa et al. 1996; Perez Garcia & Rodriguez Espinosa 2001). That the 
mid-IR emission is dominated by dust is also supported by the recent identification of silicate dust 
emission features in AGN mid-IR spectra (Weedman et al. 2005; Hao et al. 2005c; Siebenmorgen 
et al. 2005). In order to understand the spectral shapes, we have begun comprehensive modeling 
using both AGN and starburst dust radiative transfer models, the results of which will be reported 
elsewhere (Buchanan et al., in prep.). The sharpness of the break in some IRS spectra near 20 /um 
(e.g., NGG 4151; Fig. 8) is difficult to reproduce using simple dust models, however, even with 
several dust components. The IR emission is clearly produced by dust emission in the two objects 
with strong silicate absorption at 10 fiui. 

An alternative possibility is that the infrared emission may be non-thermal in origin, such 
as synchrotron emission. The broad-band photometric SEDs of nearby Seyfert galaxies provide 
evidence that the mid-IR emission of many Sy I's is predominantly non-thermal (Edelson, Malkan, 
& Rieke 1987; Ward et al. 1987; Carleton et al. 1987), although a non-thermal IR component is not 
always necessary to model AGN SEDs (Barvainis 1990). The detection of rapid variability at 10 /um 
in some nearby quasars indicates a non-thermal component of the mid-IR emission (Neugebauer &: 
Matthews 1999), and further supports that the Sy I's may also have a non-thermal IR contribution. 
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Comparison of the IR and nuclear 8.4 GHz radio flux densities suggests that the infrared emission 
is not a continuation of the radio synchrotron spectrum, as the IR fluxes are at least an order of 
magnitude too high. If the synchrotron were self-absorbed and the spectral turnover occurred at 
20 /um, consistent with the observed break in the power-law spectra, the component would have 
to be < 10^^ m in size, on the same scale as the accretion disk (Urry &: Padovani 1995). The 
typical spectral slope at frequencies below (i.e., wavelengths above 20 ^m) the turnover would 
have to be 0.8 to match the radio flux densities. This slope, and the typical slopes at frequencies 
above the turnover (wavelengths below 20 fiui) of ~-1.3, are not inconsistent with self-absorbed 
synchrotron emission, so we cannot rule this out as the origin of the power-law spectra. However, 
the energy density of such a synchrotron component would be very large (~9xl0^ ergcm"^), so 
extreme conditions would be required to conflne it. The total energy in the plasma would also be 
large, ('-^lO^'^ ergs), suggesting that such a component is unlikely. 

A third alternative for the broken and unbroken power-law spectra is free-free bremsstrahlung, 
however this would produce a much flatter slope {S^, oc u'^''^) (Barvainis 1993) than is observed in 
the spectra. 

A flnal possibility is that the IR emission is a continuation of the accretion disk spectrum, with 
a slope of, for example, 1/3 (Eq. 5.46, p. 92; Frank, King, &: Raine 1992), though more realistic 
accretion models will modify this (e.g., Litchfleld et al. 1989). This will be investigated in the 
future using X-ray - optical data in conjunction with the IR data. 

The results of our principal component analysis show that the flrst eigenvector produces more 
than 90% of the variance in the observed spectral shapes and is dominated by a starburst spectral 
shape. This suggests that the relative contribution of star-formation to the IR spectrum is an 
important factor determining the observed spectral shape. However, we emphasize that, while the 
starburst contribution dominates the observed differences in the IR spectral shapes, the starburst 
component does not dominate the total emission of each object. Most spectra are comprised of a 
power law plus a <20% contribution of eigenvector 1 (Figure 14). We found a signiflcant difference 
between the contribution of the flrst eigenvector to the spectra of the Sy I's and Sy 2's. This is 
consistent with the (on average) higher flux deflcits at 12 and 25 /xm in Sy 2's compared with Sy I's, 
suggesting that extended flux (i.e., associated with circumnuclear star formation) contributes more 
to the total mid-IR flux in Sy 2's than Sy I's. Our results conflrm the finding of Maiolino et al. 
(1995), based on ground-based photometry, that Sy 2's show more extended star formation than 
Sy I's. The sample of objects studied by Maiolino et al. (1995) were optically selected, and therefore 
our result is not due to an infrared selection bias. We have concluded that the extended emission 
is associated with star formation, however it has also been shown that extended IR emission, on 
scales of tens of parsecs, can also be associated with the AGN (e.g.. Mason et al. 2006; Bock et al. 
1998). Our data demonstrates a clear relationship between the Seyfert type and the IR SED shape 
(i.e., Sy I's tend to have power-law or broken power-law SEDs and Sy 2's have greater starburst 
contribution to their SEDs), but there is not a 1-to-l correspondence between the optical Seyfert 
type and the IR spectral shape. It is clear that star formation and complex dust obscuration play 
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key roles in producing the observed variety of spectral shapes in the IRS spectra. Homogeneous 
optical spectra, to ensure uniform Seyfert classifications, and simultaneous modeling of the optical 
and IR SEDs are necessary to elucidate these relations. 



5.2. Properties of the obscuring medium (comparison with torus models) 

Comparing the optical spectral types of the IR spectral groups (§4.2), it is not immediately 
apparent whether or not the IR spectra support orientation-dependent obscuration. Torus models 
predict silicate emission in type 1 (face-on) objects and silicate absorption features in type 2 (edge- 
on) objects (Rowan-Robinson & Crawford 1989; Pier & Krolik 1992; Siebenmorgen et al. 2004, 
though see also Nenkova et al. 2002). We would expect edge-on objects to show the most dust 
absorption. The two objects that show strong silicate absorption have optical spectral types Sy 1.9 
and HBLR Sy 2, which is not inconsistent with the unified scheme. The majority of broken power- 
law spectra that may show silicate emission features are Sy 1.2 and Sy 1.5, although the group 
includes several Sy 2's as well. The relative contribution of a starburst to the IR spectrum appears 
to be greater in Sy 2's, possibly due to a selection effect (see §5.3), so that it is difficult to discern 
any orientation-dependent difference in the AGN component of the spectra. 

The ratio of IR to radio flux densities is consistently higher for the type 1 objects in our sample 
than the type 2 objects, at all observed wavelengths. Interpreting this in the context of the unified 
scheme, this may indicate that the obscuring material is optically thick up to at least 30 /im. 
However, optically thick smooth-density torus models predict much larger differences between the 
face-on and edge-on Seyferts than the observed factor of ~2, seen in our data at wavelengths longer 
than ~15 ;um (Figure 16). For example. Pier & Krolik (1992) show that, for an optically thick torus, 
the face-on objects are expected to be at least a factor of 5 and up to several orders of magnitude 
brighter than the edge-on objects. Lutz et al. (2004) first noted that torus models cannot explain 
the observed lack of anisotropy in the IR emission of Seyferts, when they compared the 6 ^um 
emission of Sy I's and 2's, normalized to the intrinsic hard X-ray emission. Our findings, while not 
identical to those of Lutz et al. (2004) at 6 //m (discussed below), confirm that observational data 
do not match the predictions of optically thick smooth-density torus models. 

Our results are consistent with the torus becoming optically thin in the mid-IR, at around 
15 ^m. However, a mid-IR optical depth <1 would imply column density < 10^^ cm~^ (e.g., Lutz 
et al. 2004), which is inconsistent with X-ray observations of many Seyfert 2 galaxies (e.g., Bassani 
et al. 1999). The observed high X-ray column densities of Sy 2's could be reconciled with optically 
thin IR emission, if the X-ray absorption is dominated by gas in which the dust has been destroyed, 
e.g., close to the nucleus within the dust sublimation radius. There is some evidence that this may 
be possible: variability in the X-ray column densities of two Seyfert 2's on timescales of less than 
a day indicates that absorbing columns of up to 5x 10^^ cm~^ exist on scales of less than a few 
hundred Schwarzchild radii (Risaliti et al. 2005; Elvis et al. 2004). Thus our IR/radio fiux density 
ratios could indicate that the torus is optically thin above 15 ^m. As noted in §4.4 the difference in 
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the ratio of IR/radio flux densities between Sy I's and Sy 2's is particularly apparent around 10 /im; 
the 10 /um flux density from an optically thin torus is predicted to show more anisotropy than the 
surrounding continuum (Pier & Krolik 1992), suggesting that our results are also consistent with 
an optically thin torus at shorter mid-IR wavelengths than 15 /um. However, it is not clear whether 
the 6.2 fiui PAH feature associated with star formation is contributing to the apparent peak at 
10 fiui in Figure 16(b). Certainly, the steep rise in this flgure towards the shortest wavelengths 
(from 7 down to 5 ^m) suggests that the obscuring material is optically thick at these wavelengths. 

Clumpy torus models can accommodate much less anisotropy of the mid-IR emission than 
models with smooth radial density distributions (e.g., Nenkova et al. 2002). The torus emitting 
near-isotropically in the mid-IR, but with anisotropic X-ray absorption, can be explained by clumpy 
torus models, if the radial distribution of optically thick dust clumps is steep (Nenkova et al. 2006, 
in preparation). The similar ratios of IR/radio flux density of Sy I's and Sy 2's observed in our 
sample therefore may be explained by an optically thick clumpy torus. 

Our IR/radio ratios at 6 ^m appear to be in disagreement with those of Lutz et al. (2004), who 
found no significant difference between Sy I's and Sy 2's in their ratio of X-ray/IR flux densities. 
These authors used the intrinsic (i.e., absorption-corrected) X-ray flux densities as a measure of the 
AGN luminosity, and determined the ratio of the X-ray flux density to the 6 /um AGN continuum 
flux density, derived from ISO spectra. Lutz et al. (2004) find the average X-ray/IR flux ratio of 
Sy I's in their sample to be a factor of ^^2 less than the average ratio for type 2 Seyferts, which 
is contrary to the expectations of torus models in the unified scheme. Our results show a factor 
of ~7 difference at this wavelength (Fig. 16), in the opposite direction. There are several factors 
which are likely to contribute to the disagreement between our results and those of Lutz et al. 
(2004). Firstly, the relative AGN and starburst contributions to the IR fiux densities are likely 
to be different in the two studies. The ISO data of Lutz et al. (2004) were obtained using a 24" 
aperture, which will include a large contribution of extended emission; the authors account for this 
using spectral decomposition to determine the AGN continuum. In contrast our data are relatively 
small aperture and so better isolate the nuclear emission, but we have not attempted to remove 
the starburst contribution to the IR emission, which is clearly different for the type I's and 2's in 
our sample. Secondly, the sample studied by Lutz et al. (2004) was heterogeneous and includes, 
for example, Sy I's and Sy 2's with different distance and intrinsic luminosity distributions. The 
redshift and intrinsic (radio) luminosity distributions of the type I's and 2's in our homogeneous 
sample are similar (see §2) and so will not have biased the comparison of the mid-IR emission. 
Finally, the studies are using measurements at different wavelengths to represent the intrinsic 
AGN luminosities. Large corrections (up to an order of magnitude) are required to determine 
the intrinsic X-ray fiux densities from the observed fiux densities, and therefore there may be 
considerable uncertainty in the obtained intrinsic AGN luminosities, although we note that Lutz et 
al. (2004) find no evidence for a systematic offset in the absorption-corrected values which would 
bias their result. Given the various systematic biases in each of the studies, it is perhaps not 
surprising that the exact ratios of Sy 1/Sy 2/ IR emission measured are not the same. However, 
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despite this disagreement, our results at the longer wavelengths lead us to the same conclusion as 
Lutz et al. (2004): the data are not consistent with optically thick, smooth-density torus models 
(e.g, Efstathiou & Rowan- Robinson 1995), which predict strongly anisotropic mid-IR emission. 



5.3. Selection effects in the IR 

In the context of the unified scheme, the IR/radio flux density ratios of the Sy I's and Sy 2's 
in our sample suggest that the torus may be optically thick at wavelengths shorter than ~15 /nm 
(see Figure 16 and §5.2). The Seyfert sample was selected using IRAS fluxes at 12 ^m, therefore we 
might expect there to be orientation-dependent selection effects biasing our sample, that is, that 
dust obscuration will bias the sample towards type I's. A second possible selection effect, if Sy 2's 
have more star formation than Sy I's, is discussed below. To fall above the flux limit of the sample, 
the Sy 2's must contain intrinsically brighter AGNs than the Sy I's, because the Sy 2's suffer more 
MIR dust obscuration. However, we find no significant difference between the AGN luminosities 
of the Sy I's and Sy 2's in the sample (as measured by the optically thin radio emission, see §2). 
Further, if the Sy 2's have more luminous AGN emission than the Sy I's due to a selection effect, we 
may expect that the contribution of the first eigenvector would be different for high and low radio 
power sources. Again, we find no significant difference in the distributions of relative eigenvector 1 
components for high (> 10^^ WHz~^) and low power radio sources in the sample (Figure 17). It is 
possible that small number statistics are preventing us from detecting a significant difference in the 
radio powers of the Sy I's and 2's (see Table 1). The lower IR/radio flux density ratios observed 
in Sy 2's compared with Sy I's are consistent with more luminous radio sources in Sy I's. 

A second selection effect could be operating in the opposite direction to that associated with 
orientation-dependent obscuration of the nucleus. Sy 2's with intrinsically faint nuclei could make 
it into the sample if they have an additional starburst contribution to the 12 fim flux density. This 
would result in the sample being biased towards Sy 2's with a greater relative and absolute starburst 
contribution than the Sy I's in the sample. It is clear from our color-color diagram (Fig. 9) and 
principal component analysis that the observed Sy 2's show a stronger starburst contribution to 
their infrared spectra than the Sy I's, consistent with such a bias. A near- infrared study has shown 
that the Sy 2's in the 12 ^m sample are not biased towards those with powerful nuclear starbursts 
compared with optically-selected Sy 2's (Imanishi & Alonso-Herrero 2004), but the 12 ^um-selected 
Sy 2's may be biased towards more luminous extended starbursts due to the large-aperture IRAS 
observations used to select them. 

Thus we cannot exclude the possibility that the two selection effects described above are 
canceling each other out to give the same distribution of infrared luminosities for Sy I's and 2's. 
However, it is puzzling that, if the sample is biased towards Sy 2's with intrinsically brighter AGN 
and/or those with greater starburst contribution to their IR emission, that we see no differences in 
the AGN (radio) luminosities. The similar distribution of intrinsic AGN luminosities (radio powers) 
may imply that the effects of obscuration and starburst contribution are small compared with the 
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AGN luminosities. Alternatively, the statistical results may become clearer when more of our data 
are available. Further leverage on selection effects and biases can be obtained using IRAC imaging 
data, which we will address in future work. 

The selection of Seyfert samples based on the isotropic luminosity of the AGN and which 
are unbiased with respect to orientation is essential to study the physics of active galaxies and 
the AGN/starburst connection. An unbiased sample is difficult to obtain, however. The optical 
and 0.1 - 10 keV X-ray continuum fluxes suffer uncertain orientation dependent obscuration so 
favor Sy I's and discriminate against Sy 2's (e.g.. Green, Schmidt & Liebert 1987; Heckman et al. 
1995). Optical line emission, such as [O III] A5007, originates in the narrow line region so should be 
isotropic and representative of the intrinsic AGN luminosity. However it has been shown recently 
(Hao et al. 2005b) that, while the numbers of I's and 2's are comparable at low [O III] luminosity, 
Seyfert I's outnumber the 2's by a factor of about 3 at high [O III] luminosity, suggesting that [O III] 
luminosity may not be an unbiased criterion. It has further been argued, for powerful AGN, that 
some [O III] emission may come from within the torus, and therefore be orientation-dependent (di 
Serego Alighieri et al. 1997). Spinoglio & Malkan (1989) argued that selection effects are minimal in 
the 7-12 fiiQ region. The extended 12 /xm sample of Seyfert galaxies was thus selected to minimize 
wavelength-dependent selection effects. Our results have shown, however, that the nuclear obscuring 
material may be optically thick at MIR wavelengths up to 15 /um and so selection effects are still 
important in the mid-IR. The torus presumably becomes optically thin at far-IR wavelengths, 
but at these wavelengths it becomes difficult to separate nuclear emission from circumnuclear star 
formation, which emits strongly at FIR wavelengths due to the presence of cold dust. Similarly, 
at low radio frequencies and spatial resolutions, such as the FIRST survey at 1.4 GHz (White et 
al. 1997), the relatively strong contribution of star formation to the nuclear spectrum may bias 
the selection. To obtain a sample of Seyfert galaxies unbiased by possible orientation effects, it is 
therefore necessary to select galaxies based on optically thin, compact radio emission associated 
with the active nucleus (e.g., Thean et al. 2001). A large area (10 square degrees), sensitive (0.2 
mJy/beam) survey at high spatial resolution (0.2 ") at high frequency (8 GHz) would seem to be 
a useful way to select nearby Sy I's and Sy 2's unbiased by orientation effects. 



6. Conclusions 

We present the results of Spitzer IRS nuclear spectra for the first 51 of a complete sample of 
12 iiui selected Seyfert galaxies. We find the following main conclusions: 

• The spectra clearly divide into groups based on their continuum shapes and spectral features. 
The largest group (47% of the sample of 51) shows very red continuum suggestive of cool dust and 
strong emission features attributed to PAHs. Sixteen objects (31%) have a power-law continuum 
with spectral indices 05-20 um =-2.3 - -0.9 that flatten to a20-35um =-1-1 - 0.0 at --^20 ^m. Clear 
silicate emission features at 10 and 18 ^m are found in two of these objects (Mrk 6 and Mrk 335). 
A further 16% of the sample show power-law continua with unchanging slopes of 05-35 um =-1-7 
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- -1.1. Two objects are dominated by a broad silicate absorption feature. One object in the 
sample shows an unusual spectrum dominated by emission features, that is unlike any of the other 
spectra. Some features are clearly related to the starburst contribution to the IR spectrum, while 
the mechanisms producing the power-law continuum attributed to the AGN component are not yet 
clear. The precise correspondence between the IR types and the optical spectral types is unclear, 
firstly because the exact contributions of the AGN and starburst to the IR emission are not clearly 
separated by colors or principal component analysis, and secondly because the Seyfert types are 
heterogeneous and often ambiguous. 

• Principal component analysis suggests that the relative contribution of starburst-heated dust 
emission to the SED is the dominant cause of variance in the observed spectra, although all the 
spectra are dominated by a power-law component. The close detailed agreement of eigenvector 
1 with the spectrum of a starburst at A <20 /um makes it almost certain that eigenvector 1 is 
tracking the starburst component in these AGNs. Moreover, that eigenvector 1 is linked to 12 ^um 
IRAS/Spitzer flux deficit, shows that this star formation is extended and so not a purely nuclear 
phenomenon. 

• We find that the Sy 2's typically show stronger starburst contributions in their IR spectra than 
the Sy I's, confirming previous results found using photometric data. This result is in disagreement 
with the predictions of unified schemes. The stronger starburst contribution in Sy 2's compared 
with Sy I's may be a selection effect in our sample, though we find no evidence supporting the 
presence of a bias in our sample. 

• We confirm the previous finding that Sy I's show higher ratios of IR/radio emission than Sy 2's 
at ~10 fim, by a factor >6, however we find that the difference between type I's and 2's decreases 
to a factor of ~2 beyond -^15 ^m. In the context of unified schemes, this presents a challenge for 
dusty torus models, which generally predict much larger anisotropy. The observed factor of ~2 
difference between the type I's and 2's in their IR/radio ratios above 15 fim requires the standard 
smooth-density torus models to be optically thin at these wavelengths. However, the resulting low 
torus opacity requires that the high observed columns detected in X-ray absorption be produced 
in gas with very low dust to gas ratio (perhaps within the dust sublimation region). On the other 
hand, our observations may be consistent with clumpy torus models (e.g., Nenkova et al. 2002) 
containing a steep radial distribution of optically thick dense clumps. 

• The selection of our sample at 12 /im, where the torus may be optically thick, implies there may 
be orientation-dependent biases in the sample, however we see no evidence for these in our results. 
Detailed modeling of the continuum emission is underway to separate in detail the starburst and 
AGN contributions to the IR spectrum in order to place constraints on the opacity and geometry 
of the nuclear obscuring material, and to compare the relative starburst contributions of Seyfert 
types 1 and 2. 
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improvements to the paper. This work is based on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology 
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